Voter model on a directed network: Role of bidirectional opinion exchanges 
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The voter model with the node update rule is numerically investigated on a directed network. 
We start from a directed hierarchical tree, and split and rewire each incoming arc at the probability 
p. In order to discriminate the better and worse opinions, we break the Z2 symmetry {a = ±1) by 
giving a little more preference to the opinion a = 1. It is found that as p becomes larger, introducing 
more complicated pattern of information flow channels, and as the network size A^ becomes larger, 
the system eventually evolves to the state in which more voters agree on the better opinion, even 
though the voter at the top of the hierarchy keeps the worse opinion. We also find that the pure 
hierarchical tree makes opinion agreement very fast, while the final absorbing state can easily be 
influenced by voters at the higher ranks. On the other hand, although the ordering occurs much 
slower, the existence of complicated pattern of bidirectional information flow allows the system to 
agree on the better opinion. 

PACS numbers: 87.23. Ge, 89.75.Fb 



Recently, social phenomenon of the opinion formation 
has been studied through the uses of simple prototypi- 
cal models such as the Axelord's culture dissemination 
model m and the voter model j3-ll(]l|. In particular, the 
voter model on complex networks has been popularly 
studied 0-01 1 reflecting the broad research interest in 



dynamic behaviors of complex networks 11| . In the lan- 
guage of statistical physics, the voter model has the same 
symmetry as the Ising model if the Ising spin is replaced 
by the opinion of a voter who can have two different opin- 
ions cr = ±1. Likewise, the opinion formation dynamics 
in the voter model corresponds to the coarsening pro- 
cess in the Ising system, although the former is driven 
by the interface noise while the latter by the surface ten- 
sion in d-dimensional lattice [8|, resulting in interesting 
differences 0, Q. The analytic approaches have been 
successfully applied for the voter model in the regular 
d-dimensional lattice [m| and in the complex networks [^ 
13]: in d-dimensional regular lattice, the density p of ac- 
tive bonds connecting opposite opinions decays to zero 
only for d < 2, whereas for d > 2 the system never ap- 
proaches completely ordered absorbing state @, 0] ■ The 
voter model on fractal structures with d < 2 also exhibits 
the complete ordering ,10i] . The ordering dynamics of the 
voter model in various complex networks share common 
features @-[a 0|- if started from randomly distributed 
opinions, the system first approaches a quasi-stationary 
state and stays there for a life time r. As time is elapsed 
further the system approaches the absorbing state char- 
acterized by p = 0. Interestingly, the life time r of the 
disordered quasi-stationary state increases with the net- 
work size N [2-S 0] , which makes the complete ordering 
in the thermodynamic limit impossible, in accord with 
the study for regular lattices in a very high dimension- 
ality tSi]. There is now growing consensus that the or- 
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FIG. 1. (a) A directed tree network with the perfect hierar- 
chical structure is initially built with the branching ratio 3 
and the number of levels L = 3. The total number of vertices 
is given by A'^ = (3^ — l)/2. (b) Each incoming edge is visited 
one by one, and then split into two edges at the probability 
p: the existing one but with the half weight and the new one 
of the half weight from randomly chosen other vertex. Solid 
(dashed) lines denote the incoming edges of the weight unity 
(half). 



dering dynamics of the voter model depends only weakly 
on the topological details of the underlying complex net- 
works j3., _7] , and that the first and the second moments 
of the degree distribution mostly determine the ordering 
behaviors for networks with negligible degree-degree cor- 
relation [3, |j| . We emphasize that most existing studies 
of the voter and other related models have assumed that 
the networks are undirected. The most recent studies for 
the ordering dynamics of the voter model have started to 
consider the directedness of networks [IJ, [lj| . 

In the sociological and biological disciplines, directed 
networks are abundant. To name a few, the world wide 
web, the email communication network, and transcrip- 
tional regulatory networks are described better as di- 
rected networks. In such directed networks, one gener- 
ally expects that dynamical behaviors such as synchro- 
nization of phase oscillators, opinion formation in the 
voter model, and the epidemic spread of diseases, exhibit 
sharp differences from those in undirected networks. In 
Ref. |6|, both the Hodgkin-Huxley model in neuroscience 
and the sociological voter model have been studied in 
the small-world type directed network: as more edges 
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FIG. 2. The order parameter m(i) = (1/A'^)(^ ctu) averaged over 500-2000 independent runs versus the Monte-Carlo time 
step t for directed network of sizes (a) i = 4 (iV = 40), (b) i = 6 (iV = 364), (c) L = 7 (iV = 1093) with e = 0.1 and (d) L = 7 
with e = 0.2 at various values of the splitting probability p. For small sizes, the voters favor the opinion of the top cr = — 1 
at any value of p. However, at sufficiently large values of p and for larger sizes, finite fraction of voters begin to have opinions 
a = 1 against the top vertex, for a very long time, (e) Undirected network of size L = 7 {N — 1093) with e — 0.2 in which all 
edges are undirected, (f) Different version of undirected network of size L = 7 {N = 1093) with e — 0.2 in which all edges are 
undirected except the ones that connect the top voter and her lower rank direct acquaintances. 



become directed, there occur clear changes in dynamic 
behaviors, and it has been found that bidirectional infor- 
mation exchange helps consensus state to emerge in the 



voter model. Sanchez et. al. [ij] found that the directed 
links on the directed small-world networks play an im- 
portant role altering the nature of the phase transition. 
These findings show that the dircctedness of networks 
strongly affects dynamical and critical behaviors of the 
system. 

In this paper we study a sociological game model of 
voters connected through the directed network structure 
(see Fig. [T]). Our directed network originates from the 
perfect hierarchical tree structure, which has been found 
to yield the besj^synchronizability for the set of identi- 
cal oscillators 15] but not for the nonidentical oscilla- 
tors [161 . In order to build directed networks, we start 
from the tree with L layers and branching ratio 6 [L = 3, 
6 = 3 in Fig. [Ija)]. We then sequentially visit each in- 
coming edge of the weight unity and with the probability 
p reduce its weight to 1/2 and make a new incoming edge 
of the weight 1/2 from a randomly chosen vertex. It is 
to be noted that our way of splitting edges conserves 
the incoming strength for each vertex and the result- 
ing directed network contains the original tree. When 
the above procedure is repeated for all incoming edges 
in the original tree, the resulting network can possess 
loops [e.g., 2^5-^8^2in Fig. [lib)] and directed 
upward edges [e.g., 8 — ?> 2 in Fig. [Ub)]. Consequently, 
our network model embeds the hierarchical directed net- 
work, which exists abundantly in human societies: Deci- 
sion makings in corporations, governments, and religious 
organizations can be described as opinion formation on 
hierarchical structures jl7|. 



After the construction of directed networks, we per- 
form the voter model simulations on top of the network 
structures. The original voter model with two opinions 
a = ±1 for undirected networks proceeds as follows: A 
vertex v is chosen at random. One of its neighbors, say 
w, is chosen randomly among w's direct acquaintances. 
The opinion a^ of the vertex v is changed to that of 
ws. The N repetitions of the above process constitutes a 
unit time step so that each voter has chance to update its 
opinion once on average. The generalization for directed 
and weighted networks is straightforward: we only need 
to change the procedure so that w is now chosen only 
from v's incoming neighbors {w — > v) and the weight of 
this chosen incoming link is used as the probability of 
acquaintance of v and w. Suppose that we play such a 
voter model on the directed tree structure in Fig. [T][a) . 
The final absorbing state of the whole system is solely 
determined by the opinion of the top vertex since it does 
not have incoming edge and thus no one can persuade 
him otherwise [18]. In other words, if the top vertex has 
made somehow an incorrect decision, the whole society 
of voters cannot escape from the disastrous situation in 
which all voters have the same incorrect opinion of the 
top vertex. However, in a real experience, one observes 
that it is plausible that a large scale of consensus of opin- 
ions can be made although it is against the opinion of 
the top. In this regard, we give a distortion in the orig- 
inal voter model as follows: (i) the top voter is assigned 
the opinion cri = — 1 which we assume as an incorrect 
opinion (or poorly decided opinion) . All other voters are 
assigned either 1 or —1 randomly at the initial stage, (ii) 
A voter v is chosen at random, and one of its incoming 
neighbor w is also randomly picked. At the probability 



of the weight (either 1 or 1/2 in our model; see Fig. [T]) 
for the hnk connecting the two, the following procedure 
is performed, (iii) If cru,(i) ~ — 1 at time t, v changes its 
opinion to w's one, i.e., av(t+ 1) = —1 at the probability 
1 — e. On the other hand, if a^it) = 1, v always follows 
w's opinion, i.e., ay{t + 1) = 1 at the probability unity. 
In words, the better opinion cr = 1 is always accepted 
while the slightly worse opinion cr = — 1 is accepted with 
a bit of hesitation. Within this modified voter model, the 
original Z2 symmetry is broken (compare with Ref. |8|) by 
the asymmetry parameter e and we drive the system to 
be in favor of the opinion cr = 1, by giving each voter a 
slight preference toward cr = 1. The incoming strength, 
i.e., the sum of the weights of incoming arcs, is conserved 
in our model, which we interpret as that the total time 
or resources each voter can use to make contact is fixed 
to a constant. We have also used the unweighted version 
of our modified voter model in which the above step (iii) 
is performed always for the pair {v,w), regardless of the 
weight of the link. We have observed only insignificant 
differences and all the results presented in this Brief Re- 
port are for the voter model with the edge weight taken 
into account. 

The motivation of our voter model is to mimic real- 
world situations in which there are two competing op- 
tions or choices and one is intrinsically better than the 
other. The worse option can still pervade the whole sys- 
tem, if a dominant agent strongly advocates it. In our 
model, the intrinsically better option is described by the 
symmetry breaking parameter e > 0, and the strong ad- 
vocate of the worse option is represented as cr = —1 as- 
signed for the top voter. The purpose of our Brief report 
is then to reveal what can be the effects of the connec- 
tion structure to change the situation so that the better 
option cr = 1 manages to spread across the whole society. 

We construct directed networks of various sizes as de- 
scribed above and perform the simulation of our modified 
voter model. During the simulation, we measure the or- 
der parameter defined as m{t) = {^/N){J2y<^v{t)) with 
(• • •) being the average over network realizations and ini- 
tial configurations. The consensus state with completely 
agreed opinion is characterized by either m = 1 or = — 1 . 
When the splitting probability p in Fig. [1] is sufficiently 
small, we expect that the system approaches the consen- 
sus state in which all voters have opinion cr = — 1 since 
the unidirectional information fiow is abundant, and very 
few voters can disobey the opinion of the top voter. In 
contrast, as p is increased further and thus as lots of di- 
rected loops are newly formed, one expects that the situ- 
ation reverses and a large number of voters will favor the 
opinion cr = 1. In Fig. [21 the order parameter is displayed 
as a function of t for various sizes (a) L — 4 {N — 40), 
(b) L = 6 (iV = 364), and (c) L ^ 7 {N = 1093) at 
e ~ 0.1, and for comparison (d) L — 7{N — 1093) at 
e — 0.2 in directed networks, and the undirected cases 
of (e) and (f) of networks of size L = 7 {N = 1093) at 
e = 0.2 are discussed later. It is clearly seen that the 
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FIG. 3. Relaxation time r versus (a) p and [(b) and (c)] A*'. 
(a) r monotonically increases with p. (b) For e = 0.0, at 
higher values of p, r grows with A*' algebraically, whereas at 
lower values of p r grows with p logarithmically. The inset 
exhibits the result of r vs. A'^ for directed tree and undirected 
tree (i.e., where p = 0) and lines are results from the least- 
square-fit to algebraic and logarithmic functions, respectively, 
(c) For e = 0.05 and 0.1 at higher values of p, r increases with 
A^ faster than algebraic. 



above expectation becomes indeed correct as the system 
size is increased. For small sizes, the system eventually 
approaches to the absorbing state to = —1 at any value 
of the splitting probability p. However, as the size be- 
comes larger [see Fig.[2](c) and (d) for L = 7 (TV = 1093)] 
m increases toward positive values at higher values of p , 
indicating that there are more voters with cr = 1 than in 
the initial state, and m stays at an almost constant value 
for a very long time 19]. In the later stage, m appears 
to decay to —1 eventually, very slowly. Furthermore, the 
larger value of e makes more voters agree on cr = 1 at 
lower values of p as shown in Figs.[l](c) and [21(d). The 
comparison of the time decay behaviors for different sizes 
makes us conclude that for even bigger sizes, the system 



will stay at positive values of m{t) in an extremely long 
time. Due to the limitation of computing resources, it is 
difficult to decide whether m{t — ^ cx)) = 1 or — f for large 
enough value of p. 

Another interesting observation one can make from 
Fig. [2] (c) and (d) is that the approach toward m = —1 
for p = occurs much faster than toward a positive 
value of m for p — I. This indicates that although 
the system can achieve the state in which majority of 
voters keep the better opinion, it takes longer time to 
get there. We next compute the active bond density 
Pit) = E.E,evr(l-'^«'^^-)/2]/E.l^rL where vt 
(f ™) is the set of the incoming nearest vertices (edges) 
of the ith voter and |f ™| is the number of edges in the 
set £?" . The relaxation time to reach the absorbing state 
p{t) = or m{t) = — 1 is also calculated. We observe that 
as p is increased, r also increases monotonically, which 
indicates that the complete opinion agreement to ct = — 1 
takes longer time as the network structure becomes more 
different from the hierarchical structure, and r increases 
more rapidly for the higher value of e, as shown in Fig. [3] 
(a). For e = 0.0, r increases in a logarithmic way with 
N a,t p = 0.0, whereas in undirected tree the relaxation 
time increases algebraically with N as displayed in the 
inset of Fig. [3] (b), implying that the consensus of opin- 
ions occurs much faster in directed tree networks than 
in undirected ones. It is notable that the crossover from 
the logarithmic behavior to the algebraic one exists as p 
is increased in the case of e = 0.0 [see Fig. |3](b)]. Conse- 
quently, we conclude that the bidirectional information 
flows make the agreement of opinion slower through the 
formation of a large number of voters against the opinion 
of the top voter. Figure [3] (c) shows that when e > 0, r 
increases very fast with N for higher values of p (faster 
than the algebraic increase). 



For comparisons, we also study the undirected network 
version of our model in which the network is constructed 
in the same way as for directed case (see above) and then 
all edges are assumed to be symmetric, i.e., bidirectional. 
We observe that the final consensus state m{t —^ oo) > 
is reached as far as e > 0, regardless of the splitting 
probability p, which leads us to conclude that our voter 
model in directed network and undirected network ex- 
hibit a quite different behavior [see Figs. [2] (e) and [2] 
(f)]. Nevertheless, our main conclusion of the role of the 
bidirectional information flow is still valid in the undi- 
rected networks: the more abundant bidirectional opin- 
ion exchanges help the majority of voters to agree on the 
better opinion. 

In summary, we have investigated the simple sociolog- 
ical model of voters on directed networks. Starting from 
a directed tree structure, each edge weight is halved and 
a new incoming edge with a half weight is newly con- 
nected at the probability p. As more edges are split, 
making the network possess more complicated pattern of 
information flow channels, it has been found that finite 
fraction of voters can hold opinions against the opinion of 
the root voter at the top. This trend becomes more evi- 
dent as the asymmetry between the two opinions cr = ±1 
becomes bigger, making one opinion more favorable than 
the opposite, and also as p and iV become larger. The 
approach toward the absorbing state has been found to 
be very slow for large values of p. It has also been found 
that the relaxation time increases very fast with the size 
N. 
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